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ABSTRACT 

We present theoretical thermally-averaged rate constants for vibrational and rotational 
(de— )excitation of the Hl^ ion by electron impact. The constants are calculated using 
the multi-channel quantum-defect approach. The calculation includes processes that 
involve a change |AJ| ^ 2 in the rotational angular momentum J of H'^ . The rate 
constants are calculated for states with J ^ 5 for rotational transitions of the 
ground vibrational level. The thermal rates for transitions among the lowest eight 
vibrational levels are also presented, averaged over the rotational structure of the 
vibrational levels. The conditions for producing non-thermal rotational and vibrational 
distributions of in astrophysical environments are discussed. 

Key words: ISM: molecules, plasmas, molecular processes, molecular data 



1 INTRODUCTION 

Rotational and vibrational excitation of small polyatomic 
ions by electron impact is one of the important processes 
occurring in a neutral molecular plasma. In particular, the 
probability of rovibrational (de-)excitation in electron-ion 
collisions can be relatively high. The high probabilities and 
correspondingly high rate constants are driven by the non- 
Born-Oppenheimer coupling between electronic and rovibra- 
tional motions of the ion-electron system. In certain small 
1 polyatomic molecules, the coupling is particularly strong, as 
is the case for the H;^ ion. 



metastable neutral Ha molecule. Astrophysically, electron 
impact excitation of molecular ions has been observed to be 
the domina nt coUisional excitation pr ocess in some environ- 
ments fe.g. Ijimenez-Serra et al.ll2006f ). 



Owing to its importance in interstellar sp ace (|Okal 

l2006h . planetary ionospheres (^Miller et al.l ^2 000l) and lab- 
oratory exper iments (JLarsson 2000; Plas ifet al.l |2002| : 
IJohnsenI |2005| ). the Hg" ion has been studied for many 
years. In particular, processes involving electron scat- 
tering from the ion have been recently studied ex- 
perimentally and theoretically. Such processes include 
electron-impact rovibratio n al exc itation ((paurc & TennvsonI 
l2002l : iFaure et al.l l2006al. I2OO9I'). dissociat i ve recombina- 
tion jLarssonl I2OOOI: iKokoouline et al.l I2001I : Ijohnseni l2005l : 
iFonseca dos Santos et a^"2007^, electro nic excitation and 
ionis ation ([G orfinkicl fc Tennyson' 2004') and photoioniza- 
tion llBorda s" fc Hel m 1991; Stephens fc Grpcnc .19 94, 1995.; 
iMistrik et ain2000l : IKokoouline fc Greend T2004al 7 of the 



Previously, cross sections and rate constants for 
a few rotational transitions in Hj" have been calcu- 
lated jFaure fc Tennvson.20o3 . l2003l : IKokoouline fc Greend 
l2003bl : iFaure et al.l l26o6al lbl). Vibrational rate constants 
have apparently not been studied extensively for this fun- 
damental ion. Here, we present thermal rate constants for 
transitions between different rotational states of the ground 
vibrational level of HJ with low angular momentum, J ^ 5. 
We also present thermal rate constants for rotationally- 
averaged transitions between different vibrational levels. 

The next section of the article briefly discusses the 
theoretical approach used in the present calculation. A 
detailed description of the approac h is lengthy and has 
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already been published elsewhere ( Kokoouline fc Greend 
l2003bl : iFonseca dos Santos et all |2007| ). Therefore we only 
sketch here the main ideas of the approach. In section [3l 
we present the rates for vibrational (de-)excitation of Hg". 
Section |4] is devoted to the calculation of rotational rate 
constants for transitions within rotational manifold of the 
ground vibrational level. Astrophysical implications are dis- 
cussed in Section [5l Finally, Section [6] presents our conclu- 
sions. 
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2 THEORETICAL APPROACH 

The theoretical model employed in the present study is 
based on quantum defect theory. I t is discussed in de- 
tail b y Kokoouline & Greene ('2003a') , 'Kokoouline fc Greene 
1 2003b ') , Faurg_ct_al. ( 2006a ) and Fonseca dos Santos et al. 
(20071). Here, we only mention the main ideas used in the 
model. 

The energy-dependent theoretical rate constant 
Qi'<_i(i?) for a transition from the initial rovibrational level 
z to a final one i' is obtained from the corresponding matrix 
element Sf,^^'(B) of the energy-dependent scattering 
operator, S''''^''{E). The main difficulty in the theoretical 
approach is the construction of the scattering matrix 
S^!l^"{E), (the indices i and i' refer to specific rovibrational 
states of the Hj^ ion; the electron angular momentum and 
its coupling with the ion to form a total angular momentum 
eigenstate are implied as well, but these will be suppressed 
in our notation since they are diagonal quantum numbers 
in the present approximation). For our discussion, it is 
convenient to represent the index i as rv, where r = (JK) 
specifies the rotational quantum numbers, i.e. ionic angular 
momentum J and its projection K on the molecular 
axis; and v = {vi , v^2 } specifies the vibrational quantum 
numbers in a normal mode classification. We note that our 
model neglects the explicit coupling between rotational 
and vibrational angular momenta which occurs for excited 
vibrational states with h > and whi ch leads to a more 
comp licated set of quantum numbers (|Lindsav fc McCalil 
I2OOII ). 

The construction of the scattering matrix S^l^^" (E) be- 
gins from the ab initio potential surfaces of the ground elec- 
tronic state of the ion and several exc ited states of the neu- 
tral H3 molecule (|Mistrfk et al.ll2000l) . {U+{Q) for the ion 
and Un{Q) for the neutral molecule). We will use the sym- 
bol Q to specify collectively the three internuclear distances. 
For a given geometry Q, the electronic wave function of the 
outer electron of the H3 excited states resembles the elec- 
tronic wave function of the hydrogen atom. However, due to 
fact that at short distances from the ionic core the electron- 
ion interaction is different from that in the hydrogen atom, 
the electron binding energy U^{Q) — Un{Q) is generally dif- 
ferent than the corresponding binding energy l/(2n^) in the 
hydrogen atom (n is the principal quantum number). The 
departure of the f/+(Q) - U„{Q) difference from l/(2n^) 
is written as l/(2(n — /i)^), where the quantum defect is 
only weakly dependent on the principal quantum number n. 
When the energy E of the electron-|-ion system approaches 
the ionization limit (n — > 00) and becomes larger, the prin- 
cipal quantum number n looses its physical meaning, but 
the quantum defect fi does not: it gives the collisional phase 
shift 5{E) — irjj, in terms of the scattering phase in e lectron- 
ion c ollisions at energies above the ionization limit (|Seatonl 
1 19661 ). The phase shift also depends weakly on the energy E 
and it determines the scattering matrix S{E) = exp(2z(5). 
This is the reason why the energies of excited electronic 
states in the neutral molecule can be used to obtain the 
quantum defect and to describe collisions between the ion 
and the electron. 

The preceding discussion assumes that the electron 
scatters from (or is bound to) the molecular ion, which stays 
at a given geometry, i.e. the nuclei remain fixed through- 



out. This approximation is only reliable on a time scale 
much shorter than the period of ionic vibrational motion. 
In this limit, the quantum defects are functions of geome- 
try Qi m(Q)- Generally speaking, the dependence of /i(Q) is 
much smoother than the dependence of U^{Q) and U{Q). 
Although the quantum defect depends only weakly on the 
principal quantum number n, it usually depends strongly 
on the angular momentum I and on its projection A on the 
molecular axis of the ionic core. For a nonlinear triatomic 
ion such as Hg", the body- frame quantization axis is chosen 
as the normal to the plane containing the nuclei. At large 
I ^ 2, the quantum defects become small, because the elec- 
tronic wave function of high I approaches the unperturbed 
wave function in the hydrogen atom. 

Therefore, for the description of an electron-ion scat- 
tering process involving low values of /, one must obtain the 
quantum defect functions ^ia{Q) for various I and A and 
the corresponding scattering matrices S{Q) — exp {2iSiA)- 
The calculated potential energy sur faces are obtained i n the 
Born-Oppenheimer approximation l|Mistrfk et al.llioOCj '). i.e. 
in which the coupling between electronic and nuclear mo- 
tion is neglected. The quantum defects and the scattering 
matrices obtained from the potential surfaces as described 
above, therefore, fail to account for the coupling. However, 
in the H3" + e~ case a strong non-adiabatic Jahn- Teller in- 
teraction must be accounted for in order to appropriately 
describe the scattering process. It mixes the two tt elec- 
tronic states of the same orbital angular momentum I of 
the incoming electron and can be included in the scattering 
matrix b uilt from quantum defects using the formalism sug- 
gested in IStaib fc Domck^ ^9^)- 

The resulting scattering 
matrix is now not diagonal over the A quantum numbers. 
The strongest coupling is between the pvr electronic states. 
In our treatment we include only p-states (tt and a). This 
implies that thi s body-frame scattering mat rix S^i^j^{Q) is 
a 3 X 3 matrix (|Kokoouline fc Greene! [2003bl ). 

The scattering matrix constructed in this way repre- 
sents the electron-ion scattering only if the ion stays at the 
same configuration Q during the entire process, which is 
not a valid description. The physically meaningful scatter- 
ing matrix must describe the amplitude of scattering from 
a particular rovibrational state of the ion to another one, 
including the possibility of nuclear motion during the full 
collision process. We denote this as the space-fixed scatter- 
ing matrix Sf/^, where s and s' refer to initial and final 
states of the io n. The formalism of the rovibrational frame 
transforma tion (|Atabek et al.|[l973 : IJungen fc Atabeklll977l : 
lFanolll975^ allows us to use the matrix S^^{Q) to construct 
the Sfi^g matrix. In this formalism, the two matrices are con- 
sidered as two equivalent forms of the same scattering opera- 
tor in two different representation bases. The representation 
basis of 5f/^ — {s'\S\s) is the set of rovibrational energy 
eigenstates |s). The basis for the 5*^^(2) matrix is made of 
tensor products \b) = |Q)|A), where \Q) represents the vi- 
brational position eigenstates of the ion, |A) is the angular 
state vector of the electron in either the pn or pa state. 

The two representations are connected by the standard 
basis transformation formula 

sZ = J2(''\b'}{b'\S\b){b\s), (1) 

b',b 

where the summation indicates a sum over discrete indices 
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and an integration over the continuous coordinates Q and ro- 
tational coordinates (three Euler angles). Because the scat- 
tering matrix in the \b) basis is diagonal over it is con- 
venient to write {b'\S\b) as S^,'[^{Q)S{Q - Q'), the notation 
that has been already used in the above discussion. The 
explicit form of the matr ix elements (s\b) for the unit ary 
transformation is given in iKokoouline fc Greene! l|2004al lbl). 

The S^^ matrix is diagonal with respect to the rota- 
tional quantum numbers Jtot, Ktot and Mtot of the whole 
molecule and the continuous coordinate Q. Jtot, Ktot and 
Mtot are the angular momentum of the neutral molecule, 
and the two projections of the angular momentum on the 
molecular symmetry axis and on the space-fixed z-axis, re- 
spectively. To completely define the BP basis functions |&), 
in addition to A and Q, the rotational quantum numbers 
Jtot, Ktot and Aftot must be also specified; for brevity they 
are omitted in the above equation although their presence 
is implied. 

To specify the vibrational states in the SF representa- 
tion, we will use the normal mode approximation, i.e. speci- 
fying ionic vibrational eigenstates by the quantum numbers 
{vi,vi^}. The rotational part of the total wave function is 
specified by rotational quantum numbers Jtot, J, K, M, I, m, 
where J, K, and M are the angular momentum of the ion, 
and the two projections of the angular momentum on the 
molecular symmetry axis and on the space-fixed z-axis, m 
is the projection of electronic angular momentum I on the 
space-fixed 2-axis. In the following, we will not specify any 
other conserved quantum numbers that are the same in both 
bases, such as the total nuclear spin and the irreducible rep- 
resentation of the total wavefunction. 

The matrix 5*5^ obtained by the above procedure 
does not yet represent the p hysical scattering matrix 
l|Avmar et al.lll996l : ISeatonlll966h . In fact, it represents the 
actual scattering matrix S^'^^"{E) only for energies high 
enough such that all of the channels |s) are open for elec- 
tron escape, i.e. where the total energy of the system is 
higher than the energy of the highest relevant ionization 
channel \ s). When at least one channel is closed, the physical 
scattering matrix 

gphys^^^ is obtained from using the 
st andard MQDT chan nel-elimination formula (see Eg. (2.50) 



10 1(11) 



Avmar et all (Il996l ) or Eq. (38) in IKokoouline fc Greeii3 



l|2004ah ') 



In terms of the energy-dependent scattering matrix 
S^'^'"" [E), the cross-section for rovibrational (de-)excitation 
of the ion from the initial state |s) is written (in atomic 
units, a.u.) as: 



2Eel{2J+l) 



^ (2Jtot + l)|5: 



(Jtot 5^^t0t) I 2 



(2) 

where E^i is the relative kinetic energy of the ion and the 
electron before the collision. In the above expression it has 
been assumed that the initial js) and final |s') states are 
different. 



3 RATE CONSTANTS FOR VIBRATIONAL 
(DE-)EXCITATION 

If one is not interested in the rotational structure of iiu- 
tial and final vibrational states, the final cross-section (or 




total energy , with respect to { 00 }(00) (cm ) 

Figure 1. (Color online) Probabilities of rovibrational excitation 
of the H J ion calculated using the full rovibrational frame trans- 
formation. Only transitions from the ground vibrational level 
{00°} are shown. The {OO"} {01^} probabilities oscillate a 
lot below 3000 cm~^ and become less energy-dependent above. 
The oscillations are due to the strong rotational coupling between 
individual rotational levels of the initial and final states of the 
ion. When averaged over the initial and summed over the final 
rotational states and averaged over the appropriate energy dis- 
tribution, the resulting probabilities are similar in magnitude to 
the probabilities shown in Fig. [2] The labels on top of the figure 
indicate different rovibrational ionization limits. Notice that the 
zero of energy in the figure is set to the energy of the forbidden 
rovibrational level {000}(00). 



thermal rate constant) has to be averaged over the ini- 
tial rotational levels and summed up over the final rota- 
tional levels. This can be done using the full rovibrational 
frame transformation technique described above. However, 
the cross-section averaged over initial rotational levels and 
summed up over the final rotational levels is very similar 
to the one obtained neglecting the rotational structure of 
the ion. The cross-section obtained accounting for the rota- 
tional structure has more resonances due to interaction be- 
tween rotational states but the averaged value is close to the 
value obtained without the rotational structure. The quan- 
tity of interest in astrophysical applications is the thermally- 
averaged rate constant. Because the thermally-averaged rate 
is not sensitive to the position of individual resonances in the 
energy-dependent cross-section, the calculations with and 
without the rotational structure give the same result. 

Figs.[T]and[2]compare our calculations with and without 
rotational structure included, for the probabilities of rovi- 
brational and vibrational (de-)excitation of the ion. Fig. [T] 
shows in detail the rovibrational transitions from the ground 
to the first excited vibrational level {01^} with Jtot = 2 in 
para-Hg". In order to compare with the results on Fig. [2j 
one would need to take a sum over final quanta and average 
over initial J and K and account for all possible Jtot and 
Ktot similar as it is done in Eq. (2] This would mean that 
to achieve a converged result at reasonably high energy (~ 
2000 cm"'^) calculations would be needed for all Jtot up to 
10. This would require a tremendous numerical effort if the 
fully quantum approach were to be applied. 
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^000 3000 4000 5000 6000 7000 8000 9000 10000 
electron scattering energy (cm ') 

Figure 2. (Color online) Probabilities of vibrational excitation 
from the ground vibrational level {OO"} to several excited vibra- 
tional levels calculated using the vibrational frame transformation 
only. Energies of vibrational thresholds are labeled with arrows 
and the corresponding vibrational quantum numbers. 
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100 1000 
temperature (K) 

Figure 3. (Color online) Thermally-averaged rate constants for 
several (de-)excitation transitions obtained by direct integration 
using Eq. ^ (solid lines) and the approximate formula of Eq. 
10 (dashed line). The averaged probabilities for vibrational (de- 
)excitations are listed in Table [l] 



The thermally averaged rate constant ath(T) (in a.u.) 
is obtained from the energy-dependent cross-section (j[E) as 



ath{T) 



8tt 



(27rfcT)»/2 



a{E^i)e "T E^idE^i 



(4) 



where T is the temperature. Temperature dependencies 
ctth{T) for different (ro-)vibrational transitions w — >■ u' ob- 
tained using Eq. (|4]) are shown in Fig. [3] as solid lines. 

For further discussion, it is convenient to represent the 
cross-section a{Eei) in the form 



fc2 



P{Ee. 



(5) 



where k is the wave vector of the incident electron, P[Eei) 
is the probability for vibrational (de-)excitation at collision 
energy Ef-i - Figures[T]and[2]suggest that on average (here we 
mean a running average taken over a few intervals between 
resonances), the probability behaves approximately as a step 
function 



{P{E,i))=Po e{E,i - A,, , 



(6) 



where A„/„ — £„/ — E^ is the threshold energy for (ro- 
) vibrational excitation (if Ey' — E^ > 0); A„/„ = for de- 
excitation (if Eyi — Ey < Q), 6 is the Heaviside function, Po 
is a constant. The above approximation for the probability 
is accurate enough to calculate the thermally-averaged rate 
constant that is not sensitive to the detailed resonance struc- 
ture of the energy-dependent (de-)excitation cross-section. 
Using Eq. (|6} the thermally rate constant of Eq. (|4]) be- 



(7) 



The above formula with only one parameter Po provides a 
very good approximation for the actual thermally averaged 
rate constant. It is demonstrated in Fig. [3] that compares 
the thermal rate constants for different (de-)excitation tran- 
sitions obtained with the direct numerical integration using 
Eq. (|4]) and with the approximate formula of Eq. 0. There- 
fore, for practical applications, it is convenient to provide 
just averaged probabilities Po and the energies of vibrational 
thresholds Ey for each pair of vibrational (de-)excitations. 
These parameters are listed in Table [T] for all combinations 
of the first eight vibrational states of H^. Note that the 
conversion factor from a.u. to cm^/s is 6.126x10^^. 



Neglect of the rotational structure of the initial and final 
vibrational state simplifies considerably the numerical cal- 
culation. The complete rovibrational frame transformation 
of Eq. ([1} is reduced to the vibrational frame transforma- 
tion if the rotation is neglected, i.e it is carried out using the 
following formula 

SfJlv')(vA) = (^v'\S'k^,iv{Q)\'v^ , (3) 

where the brackets imply an integration over the vibrational 
coordinates only. Many elements among <S'(^f^/)(„j\) are zero 
because of the symmetry of vibrational wave functions and 
matrix elements S^,^{Q). The vibrational (de-)excitation 
II — ^ «' cross-section obtained from the scattering matrix of 
Eq. ((3]) should be averaged over A and summed over A'. 



4 RATE CONSTANTS FOR ROTATIONAL 
(DE-)EXCITATION 

If the temperature T of the Kg +e~ plasma is not very high, 
such that only the ground vibrational level {OO"} of HjJ' is 
significantly populated, knowledge of rate constants for tran- 
sitions r — > r' between individual rotational levels r and r' 
of {OO"} may be important for the analysis of experimental 
or astronomical spectra. For this purpose, we have made a 
detailed analysis of transitions between individual rotational 
states of the ground vibrational level. The calculation of ro- 
tational (de-)excitation rate constants was carried out using 
the cross-section of Eq. ^ and numerical integration of Eq. 
Q. Examples of the thermally-averaged rate constants for 
the rotational (de-)excitation are shown in Fig. |4l 
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Table 1. Parameters Pq for several vibrational transitions that can be used in the approximate formula, Eq. (0, for 
thermally-averaged rate constants. Initial states, v are given in the upper row, final states v' - in the left column. The 
upper row specifies also the energies Ey (in cm-l) of vibrational levels. The probabilities Pq are obtained by fitting 
numerical dependencies obtained by a direct integration of Eq. Q. That is why Po{v' —> v) is not exactly equal to 
Po{v — >■ v'). Notice that multiplicity factors of doubly degenerate vibrational states E are taken into account in the 
probabilities. For example, Po{v — > v') ^ 2Po{v' — > v), if the vibrational states v and v' are the states of the Ai and E 
irreducible representations correspondingly. 
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Table 2. Parameters ao, ai, a2, as of the fit polynomial Pm {x) of Eq. ^ for several transitions 
between rotational states of the E" irreducible representation of the coordinate part of the ion- 
electron system. The vibrational level of the ion is the same in the initial and final states. The 
upper line specifies the pairs (JiKi) — (J2K2) of rotational states for which the parameters are 
fit. For convenience, we also specify (second line of the table) the threshold energy A^/^ for each 
transition. All but the first rows in the table have two numbers: The upper number in each cell 
corresponds to the transition (JiKi) — > {J2K2), the lower number corresponds to the reverse 
transition {JiKi) <— (J2-R'2). 
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Figure 4. (Color online) Thermally-averaged rate constants for 
several rotational (de-)excitation transitions (JK) {J'K') of 
the Hg" ion (solid lines). The vibrational level, {OO"}, is the same 
in the initial and final state of the ion. The dotted lines show a 
few examples of the numerical fit using Eq. (O. 
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Figure 5. (Color online) P^s = 

VT cxp {A^,^^/T) a[J'K' <- J K]. Note that P^g is only 
weakly dependent on temperature. It is used to obtain a 
cubic polynomial fit Pm (ln(r)) in Eq. JO} for each transition 
(JK) {J'K') with four parameters that are listed in Tables (2] 
EIH and[5] 
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490 


1092 


602 


614 


0.0 


1.39e-5 
9.38e-6 


3.69e-6 
1.87e-6 


8.36e-6 
5.99e-6 


5.94e-6 
3.48e-6 


9.79e-6 
7.42e-6 


1.07e5 
8.01e-6 


0.1 


-2.73e-6 
-1.66e-6 


5.97C-8 
1.37C-7 


-1.34e-6 
-7.52e-7 


-5.51C-7 
-1.85e-7 


-2.51e-6 
-1.73e-6 


-7.78e-8 
3.53e-7 


02 


5.98C-7 
3.73C-7 


-2.38e-8 
-3.22e-8 


4.05e-7 
2.63e-7 


8.42e-8 
2.32e-8 


4.69e-7 
3.25e-7 


-1.77e-7 
-2.22C-7 


03 


-4.35e-8 
-2.79C-8 


1.53e-9 
1.98e-9 


-3.55e-8 
-2.46e-8 


-4.e-9 
-7.4e-10 


-3.06e-8 
-2.16e-8 


1.6e-8 
1.77e-8 



Table 3. Parameters for several transitions between rotational states of the E' 
irreducible representation of the ionic wave function. See Table [2] for details. 



As is evident from the figure, the rotational rate con- 
stants behave approximately according to Eq. 0, where 
A„/„ should be replaced with the rotational threshold en- 
ergy, A^ir- = {Eri — Er)6{E^i — Er) ■ Howcver, there is 
a weak departure from the dependence of Eq. 0. It is 
clearer in Fig. [S] where we plotted an "effective" value 
PeS = a^i ^r{T)VT exp ^— fr^^ of the parameter Pq as a 
function of ln(r). Notice that on average, the quantities 
Poft [{J'K') ^ (JK)] and P<,h[{JK) ^ {J'K')] for the two 
opposite processes are related by the principle of detailed 
balance 

(2J-|-l)Peff [{J'K') ^ {JK)] ^ {2j'+l)P,s [{JK) ^ {J'K')] , 

(8) 

because Pch{x) represents the thermally averaged probabil- 
ity of the rotational transition per one electron-ion collision 
(see also Eq. ((2|). In order to simplify eventual applications 
of the calculated numerical constants Q:r',r(r), we fitted the 
numerical rate constants representing the "effective" value 
of Po in Eq. ((Tj to a cubic polynomial of ln(T). For this 
we used the following analytical interpolation formula for 

ar',r{T) 

ar'^r{T) = exp ^--^T^^ Pm{x) , where 

Pm{x) = ao + a\x -I- aix^ -f 033;^ and (9) 
X = ln(T) , 

where the constants {i = 0, 1, 2, 3) are obtained for each 
individual transition r r' from a numerical fit. The quan- 
tity Pm{x) in the above equation has a meaning of (de- )ex- 
citation probability that varies weakly with energy. We used 
the subscript m to distinguish it from Poff (that is energy- 
independent) and to stress that Pm{x) is given by a model 
polynomial. The constants obtained, a^, are listed in Tables 
[2I [3l U and [5] for several combinations of initial r = {J,K) 
and final r' = {J'K') rotational states of Hjj'. The numerical 
values of ai listed in the table are such that, when plugged 
into Eq. they give rate constants in units of cm^/s. Tem- 
peratures in the calculation of x = ln(r) should be in K. 
Notice that (2Ji -I- l)ao[l -s- 2] (2J2 + l)ao[l 2] due to 
the principle of detailed balance. 



5 ASTROPHYSICAL IMPLICATIONS 

The vibrational and rotational distribution of ions in in- 
terstellar space, planetary ionospheres and in laboratory is 





(10)-(30) 


(30)-(50) 


A,,,(K) 


619 


1085 


ao 


9.54e-6 


6.01e-6 




3.85e-6 


3.58e-6 


ai 


-3.87e-7 


4.14e-8 




-3.32e-8 


1.64e-7 


0.1 


1.43e-7 


1.04e-7 




3.71e-8 


4.12e-8 


13 


-1.155e-8 


-1.06e-8 




-3.50e-9 


-5.23e-9 



Table 4. The table gives parameters a; for several transitions 
between rotational states of the A'^ irreducible representation of 
the ionic wave function. See Table [2] for details. 





(33)-(43) 


(33)-(53) 


(43)-(53) 


A.v(K) 


494 


1101 


607 


00 


1.39e-5 


2.16e-6 


1.26e-5 




9.92e-6 


1.31e-6 


9.84e-6 


Ol 


-1.34e-6 


7.37e-8 


-2.07e-6 




-5.56e-7 


8.18e-8 


-1.41e-6 


02 


4.56e-8 


-1.81e-8 


3.86e-7 




-5.41e-8 


-1.79e-8 


2.65e-7 


03 


3.38e-9 


1.09e-9 


-2.75e-8 




7.95e-9 


1.07e-9 


-1.95e-8 



Table 5. The table gives parameters a; for several transitions 
between rotational states of the A'^ irreducible representation of 
the ionic wave function. See Table [2] for details. 



determined by the competition between radiative and coUi- 
sional processes. While the ionization level is generally low 
in astrophysical plasmas, electrons can still play a role in 
the molecular excitation because electron-impact rates ex- 
ceed those for excitation by neutrals (H, He, H2) by several 
orders of magnitude. Detailed excitation models, including 
coUisional data for all relevant colliders, are therefore re- 
quired to derive reliable column densities from astronomical 
observations. 

An important concept in this context is the so-called 
critical density, Ucr , which is defined as the density at which 
the coUisional rate is equal to the spontaneous radiative rate. 
The usual definition refers to a specific transition in a two- 
level approach. For a multi-level system, neglecting opacity 
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T 


{01^} 


{10"} 


{02"} 


{02^} 


{lin 


{20"} 


{03^} 


10 
100 
1000 


1.3(9) 
4.0(9) 
1.1(10) 


2.4(7) 
7.7(7) 
2.0(8) 


3.1(9) 
9.7(9) 
2.3(10) 


2.6(9) 
8.1(9) 
2.1(10) 


9.5(8) 
3.0(9) 
8.9(9) 


4.8(7) 
1.5(8) 
4.8(8) 


1.3(9) 
4.2(9) 
1.3(10) 



Table 6. Critical electron density, ra^r in cm ^, at 10, 100 and 1000 K, for 
vibrational levels of Hj}'. Powers of 10 are given in parentheses. 



effects, a practical definition is to refer to a specific level a 



ncr{s,T) = 



(10) 



E,,a(s^s',r)' 

where A{s — >■ s') are the Einstein coefficients for sponta- 
neous emission, a{s — >■ s' ,T) are the coUisional rates and 
the sums run over all possible transitions s — >■ s'. Consid- 
ering electron collisions only, these latter will maintain the 
level s in local thermodynamic equilibrium (LTE) for elec- 
tron densities ^ ncr{s, T), while deviations from LTE in- 
cluding population inversions are expected for densities Ue < 
Ucris, T). For Tie <^ Ticris, T) , clectron collisions will be neg- 
ligible. Non-LTE effects caused by H2 c ollisions have been 
investigated bo th in interstellar clouds (lOka &: Epd |2004|: 
Oka et al.|[2005l ) and in the Jovian atmosphere (|Melin et al.l 



20051 ). We note in this context that microcanonical sta- 



tistical calculations have been performed recently to esti- 
mate thermal state-to-state rate coeffic ients for the Hj"-|-H2 
reaction and its deuterated variants jPark fc Lightl [20071 : 
iHugo et al.ll2009l ). Electron-imp a ct rota tional excitation has 
been considered bv lFaure et al.l (|2006bl ) but, to the best of 
our knowledge, electron-impact vibrational excitation has 
been so far ignored in non-LTE modelling. 

Eq. (|10|) was computed with the vibrational and ro- 
tational rates presented in the previous sect i ons. Ei nstein 
A coefficients were taken from Dinelli et al. ('l992al,'^ for 
vibrational transitions and from Pan & Oka (1986) for ro- 
tational transitions. Results are presented in Tables [6] and 
[T] It can be noticed that critical densities for vibrational 
and rotational levels differ by typically 6 orders of magni- 
tude: they range between 10^ and 10^°cm~^ for the for- 
mer and between 10~^ and 10*cm~^ for the latter. In 
the diffuse interstellar medium, the e lectron density is 
~0.1 cm~'^ (jSlack fc van Dishoeck|[l99ll ) while it can reach 
abo ut 10^ cm~^ in pla netary atmospheres (see for exam- 
ple [LystruE^^L^al] ^200^) . As a result, non-LTE rotational 
populations are expected in interstellar clouds whereas ro- 
tational levels should be at or close to LTE in the jovian at- 
mosphere. A non-thermal rotational distribution of was 
actually observed towards galactic center clouds where the 
metastable (3, 3) level has a populatio n comparable to that 
in (1, 1) despite being 361.5 K higher (|Oka et al.ll2005l ). On 
the other hand, electrons are expected to be negligible in vi- 
brationally exciting H;^ in the interstellar medium but they 
could establish a non-LTE vibrational population of H;^ in 
planeta ry environments , as observed in the jovian thermo- 
sphere (|Kim et al.ll 19921 ). 

Table[6]shows that the two levels {10°} and {20*^} have 
critical densities significantly lower than the other levels. 
This dir ectly refiects the lo wer Einstein A coefficients of 
~1 s~ (IDinelli et al ] |l992al '). We note that the values are 
however not low enough to ensure LTE vibrational popula- 
tion in planetary atmospheres where Ue < 10^ cm"''. 

Table [3 shows that at 10 K the levels (22) and (44) 



have much higher critical densities than the other levels. 
This actually reflects the fact that these two levels can de- 
populate coUisionally through excitation only since rates for 
rotational transitions with /S.K 7^ are null within our 
treatment. These rates, as those with |AJ | > 2 , were actu- 
ally estimated bv lFaure fc TennvsonI (|2003l ) and lFaure et al.l 
(20063) and were found to be 3 to 4 orders of magnitude 
smaller than those with AJ = ±1,±2 and l\K = 0. The 
(44) level is particularly interesting because an astrop hysical 
mase r is predicted in the (44)— )-(31) transition of (|Blackl 
I2OOOI ). It should be noted that the selection rules for the 
(forbidden) rot ational radiative transitions are AJ = 0,± 1 
and AK = ±3 (|Pan fc Okalll986l : iMiller fc Tennvsonlll988l '). 
Critical densities suggest that electrons might contribute, in 
some environments, to create and maintain the necessary 
population inversion. Note that in Table [7] the (33) level is 
not listed since it is metastable. In this case, the concept 
of critical density is meaningless. Finally we emphasize that 
critical densities provide guidance at the order of magnitude 
level and that a detailed non-LTE modelling, including all 
relevant colliders, is necessary to properly quantify devia- 
tions from LTE. 



6 CONCLUSION 

In this study we have performed calculations of thermally 
averaged rate constants for rotational and vibrational tran- 
sitions in caused by an electron impact. The calculations 
were made from the flrst principles using the quantum de- 
fect approach. The rotational rate constants are calculated 
for the ground vibrational level of the ion in the initial and 
final states. The rate constants for transitions between differ- 
ent vibrational levels are calculated neglecting the rotational 
substructure of each vibrational level, which corresponds to 
averaging over initial rotational states and summing over the 
final rotational states. The obtained thermally-averaged rate 
constants are well described by the analytical formula of Eq. 
(O with the parameter Pq, that can be considered as tem- 
perature independent for vibrational transitions, and weakly 
dependent on temperature for the rotational transitions. For 
the rotational transitions, we have made a numerical fit of 
the parameter by a cubic polynomial. The numerical values 
of the fitting procedure are provided in Tables 1-5. The pre- 
sented thermally-averaged rate constants can be useful in 
interpretation of hydrogen-dominated plasma experiments 
as well as for modelling interstellar clouds and planetary at- 
mospheres, where the H3" ion is present. The computation of 
critical densities suggests in particular that electrons could 
establish non-LTE rotational populations of Yi^ in diffuse 
interstellar clouds and non-LTE vibrational populations in 
planetary atmospheres. 

This work has been supported by the National Science 
Foundation under Grant No. PHY-0855622, by the Depart- 



8 V. Kokoouline, A. Faure, J. Tennyson, and C. H. Greene 



T 


(21) 


(31) 


(41) 


(51) 


(30) 


(50) 


(43) 


(53) 


(22) 


(32) 


(42) 


(52) 


(44) 


(54) 


10 


2.5{-l) 


2.7(1) 


4.6(2) 


2.4(3) 


3.6(1) 


2.2(3) 


3.6(1) 


2.2(3) 


3.9(15) 


1.3(1) 


1.2(2) 


1.1(3) 


6.4(23) 


1.3(1) 


100 


9.6(-l) 


9.1(1) 


1.2(.3) 


7.9(3) 


1.1(2) 


6.0(3) 


6.7(1) 


2.0(1) 


3.6(1) 


4.3(1) 


3.8(2) 


3.8(3) 


2.3(0) 


1.8(1) 


1000 


1.9^0) 


l.r>(2) 


2..5(3) 


2.7(1) 


2.1(2) 


1.8(1) 


1.1(2) 


7.0(1) 


3.5(0) 


7.7(1) 


9.1(2) 


1.3(1) 


3.5(-2) 


1.8(2) 



Table 7. Critical electron density, nor in cm ^, at 10, ICQ and 1000 K, for several rotational levels of Hg" (the levels are grouped a<;cording to 
the corresponding irreducible representations). Powers of 10 are given in parentheses. 
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